RNA Biology

Due to its crucial role in proteobacterial RNA processing and degradation, the synthesis of the endonuclease RNase E is tightly controlled by a post-transcriptional autoregulatory mechanism governed by the 5' untranslated region of rne mRNA, the transcript that encodes it. Previous studies have identified one structural component of the rne 5' untranslated region (hp2) that helps to mediate this process in Escherichia coli by binding to RNase E. Here, we show that another component (hp3), a branched stem-loop, makes an even greater contribution to feedback regulation of RNase E synthesis and that the effects of hp2 and hp3 are additive, suggesting independent function. Despite the structural complexity of hp3, a single constituent branch appears to be the principal source of its regulatory activity. This branch also enables hp3 to bind RNase E with high affinity. The secondary structure of this branch, including the presence of two internal loops essential for its impact, is conserved in diverse Enterobacterales species, whereas the sequence of the internal loops is not. We propose that hp2 and hp3 mediate autoregulation by recruiting RNase E to the rne transcript.

Transfer RNA (tRNA) is highly modified and difficult to profile using reverse-transcription-based sequencing. Here, we established a direct RNA sequencing (tRNA-DRS) workflow to simultaneously characterize tRNA abundance and modification signals in mouse liver. The workflow included small-RNA enrichment, deacylation, splint-adapter ligation, barcoding, reverse transcription, and nanopore sequencing on PromethION with modification-aware basecalling. From 1,678,666 raw reads, 1,421,279 (84.67%) reads were retained after 20-200 nt length filtering, and 17,094 uniquely mapped reads (MAPQ [≥] 10) were obtained. These reads covered 244 tRNA genes and represented all 23 tRNA isotypes. tRNA abundance was strongly skewed, with Asp, Gly, Ser, and Lys accounting for 56.6% of uniquely mapped reads. Read-body coverage indicated that most uniquely mapped reads spanned substantial portions of tRNA molecules (mean 71.1%, median 75.9%; 83.3% reads covering [≥]50% of the reference). Modification analysis identified 41,013 candidate sites, of which 711 passed stringent filtering (modification fraction [≥]50%, coverage [≥]5), spanning eight modification types and dominated by pseudouridine and m5C. Filtered sites showed position-specific clustering and isotype-dependent patterns. Together, these results support tRNA-DRS as a practical single-assay framework for simultaneous tRNA expression and epitranscriptomic profiling in mammalian tissue.

Double-stranded RNA (dsRNA) biopesticides that silence essential fungal genes by spray-induced gene silencing (SIGS) are advancing towards registration, underpinned by the claim that silencing is confined to sequence-matched target organisms, a claim that has never been tested systematically against beneficial fungi. We screened twelve dsRNA constructs (the whole transcripts of eleven canonical antifungal target genes of Fusarium graminearum: three CYP51 sterol 14-demethylases, six chitin synthases and two {beta}-tubulins, plus a concatenation of the three CYP51 transcripts modelling the flagship CYP3RNA design) for perfect 21-nucleotide (21-nt) identity, the primary match criterion of recent regulatory bioinformatics frameworks, against the transcriptomes of seven non-target fungi, including commercial biocontrol agents (Trichoderma harzianum, T. virens, Beauveria bassiana, Metarhizium anisopliae, M. brunneum), the arbuscular mycorrhizal fungus Rhizophagus irregularis and Saccharomyces cerevisiae. Eleven of the twelve constructs carried off-target hazard; only CYP51C was completely clean. TUB2b carried 580 perfect 21-nt matches and contained no hazard-free window of 250 consecutive 21-mer start positions, whereas 26.6-60.0% of screened windows were hazard-free in the designable genes. Hits landed overwhelmingly on orthologues of the targeted gene and concentrated in the Hypocreales; R. irregularis and S. cerevisiae were nearly clean. These are hazard findings, not risk findings: perfect siRNA-length matches to beneficial fungi are common in canonical antifungal dsRNA targets, and current regulatory screening panels contain no fungus that would detect them.

MicroRNAs have been likened to the "dark matter" of eukaryotic genomes, reflecting their pervasive regulatory influence. MicroRNAs were first identified through genetic studies of developmental timing in the nematode Caenorhabditis elegans. Let-7 was the first microRNA recognized to be broadly conserved. The principal target of Let-7 in the developmental timing pathway is the TRIM-NHL RNA-binding protein LIN-41. During the L4 larval stage, Let-7 represses lin-41 translation by binding to two Let-7 complementary sites in the lin-41 3'UTR. Despite the importance of microRNA-based translational regulation, the underlying molecular mechanisms are incompletely understood. Through genetic analysis, we discovered an unrecognized feature of the mechanism by which Let-7 controls lin-41 translation. This mechanism requires a 5'-regulatory exon containing a seven-amino acid upstream open reading frame (uORF) and conserved sequence elements. Genome editing indicates that the specific uORF amino acid sequence itself is not important. Our data suggest that uORF translation and 5'UTR structure limit initiation at the downstream lin-41 start codon, enabling tight control by Let-7. Without this mechanism, the Let-7 microRNA is unable to properly regulate lin-41 to enable proper development.

Multiple lines of evidence suggest that RNase H1-sensitive RNA:DNA hybrids (thereafter referred to as hybrids) could disrupt the stability and restart of stalled replication forks. The prolonged over-expression of RNase H1 is the most popular method to manipulate hybrids in situ but it cannot discriminate between direct and indirect effects resulting from their loss. To overcome this limitation, we developed iCoRD (in Cellulo RNase H Delivery) to rapidly deliver comparable amounts of either Escherichia coli RNase H1 (RnhA) or control proteins in live cells, with a negligible impact on the transcriptome. iCoRD-delivered RnhA recognises hybrids in situ, is enriched in the vicinity of active replication forks and rescues replication fork progression in different stress conditions. While iCoRD-delivered RnhA had a clear impact on stress-induced, fork-proximal hybrids and R-ChIP signals, it had no effect on DRIP or S9.6 Cut&Tag signals, even when these were induced by stress. Our results strongly suggest the existence of distinct populations of hybrids that are best mapped by different methods and are differently sensitive to RNase H1 activity in situ.

Double-stranded RNAs are inflammatory byproducts of in vitro transcription of single-stranded RNA. Here we investigate removal of dsRNA byproducts using dsRNA affinity chromatography (dsRNA-AC) and two common mRNA molecules: 1) GFP, and 2) high dsRNA wild-type firefly luciferase (fLucWT). For GFP mRNA, dsRNA levels decreased by 2,000-fold, and mRNA recovery was high. In contrast, dsRNA levels decreased only two-fold for fLucWT, and mRNA yields were lower. Biophysical characterization revealed that fLucWT contains polydisperse higher-order structures that interfered with dsRNA-AC. Sedimentation-velocity analytical ultracentrifugation experiments demonstrated the higher-order structures were driven by mRNA concentration. Decreasing fLucWT concentrations increased the effectiveness of dsRNA-AC by reducing polydispersity. dsRNA levels were measured with ELISA, dot-blots, immuno-northern blots, and a cell-based interferon release assay, enabling comparison and cross-validation of these analytical methods. dsRNA-AC is a new chromatography modality, and this study identifies sample polydispersity as a key determinant of efficient dsRNA removal.

RNA interference (RNAi) offers a sequence-specific approach to pest control. In insects, Dicer-2 processes double-stranded RNA (dsRNA) into small interfering RNA (siRNA) duplexes, from which the RNA-induced silencing complex (RISC) retains a guide strand. Only antisense-loaded RISC can mediate cleavage of the target transcript. However, how sequence features shape the RISC-bound siRNA pool in pests remains poorly understood, limiting opportunities for sequence optimization. Here, we profiled RISC-bound siRNAs following injection of 34 insecticidal dsRNAs targeting 11 essential genes in Tribolium castaneum larvae. We computationally reconstructed 7,879 siRNA pairs and examined associations between sequence features and strand bias. Differences in GC identity at terminal paired positions 1 to 5, used as a proxy for local thermodynamic asymmetry, correlated with strand bias, with the strongest correlations at the first two paired positions. ORF targeting and reduced predicted antisense self-folding were also associated with higher antisense fractions. Analysis of non-templated terminal additions revealed predominantly 3-prime uridylation, a known signature of small RNA turnover, along with putative 3-prime trimming. Among ORF-associated siRNA pairs, sense strands showed higher relative U-tailing abundance, based on 3-prime uridylated and putatively trimmed-and-3-prime-uridylated reads relative to perfect 21-nt reads, than antisense strands. Antisense strands with the least predicted self-folding also showed low relative U-tailing abundance. These observations are consistent with sequence-dependent contributions from both guide-strand selection and differential post-RISC-loading siRNA retention, although a causal link remains to be established. The identified associations provide a basis for testing whether dsRNA sequence optimization can improve pest control efficacy and reduce off-target activity.

Translation initiation is a critical stage of protein synthesis in prokaryotes, where ribosomes assemble on mRNA to begin decoding genetic information. The interactions between ribosomal RNA (rRNA) and mRNA significantly influence this process, with features like ribosomal binding site (RBS) affecting ribosomal dynamics. However, direct experimental validation of these types of interaction is yet to be provided. In this study, we developed an optimized small sub-unit sequencing (SSU-seq) protocol for Escherichia coli, enabling the detailed examination of ribosomal complexes and subunits dynamics. Unlike Ribo-seq, which primarily maps fully assembled ribosomes, SSU-seq captures both ribosomal subunits and intact ribosomes, offering a comprehensive view of translation initiation and elongation events. By integrating experimental data with predictive models, our study advances the understanding of translation regulation in prokaryotes and underscores the utility of SSU-seq in decoding the complexities of ribosome-mRNA interactions. Our analysis reveals distinct patterns of small subunit distribution across the entire mRNAs with higher density near positions with stronger rRNA-mRNA hybridization. The results may suggest that rRNA-mRNA hybridization positions along the entire transcript facilitate efficient start codon recognition by maintaining proximity between the small subunit and the mRNA.

Background: Long non protein coding (lnc)RNAs are ubiquitous transcripts known to be major regulators of cellular processes in eukaryotes. In diatoms, specific lncRNA expression is tightly correlated to nutrient fluctuations, which suggests function, but their biological relevance remains to be established. In order to fill this knowledge gap, we sought to design a protocol to generate diatom lncRNA knockout mutants by CRISPR-Cas9 enabling subsequent functional studies to shed light into their mechanistic roles. Since lncRNAs do not contain open reading frames that can be targeted by point mutations for the generation of functional knockouts, we designed a multiplexed CRISPR-Cas9 plasmid harboring two independent gRNA for the specific ablation of lncRNA loci in P. tricornutum. As proof of concept, we used our plasmid-based system to target three phosphate depletion responsive lncRNAs. Results: Using the golden gate based uLoop system, we assembled a single plasmid containing all the necessary components to generate CRISPR-Cas9 large genomic ablations of lncRNAs in P. tricornutum using bacterial conjugation. Three intergenic lncRNAs (Linc91, Linc77 and Linc24) responsive to phosphate fluctuations were selected and independently targeted. Around 32% of the exconjugants presented a monoallelic ablation of the targeted loci and 11% presented a biallelic ablation on the first PCR-screening round. Mutated lines, harboring biallelic lncRNAs ablations for each lncRNA targeted, were successfully generated for each of the targeted lncRNAs. Additionally, a first line of phenotypic analysis revealed significant growth rate differences for one of the mutants ({Delta}Lnc24) under both full and phosphate depleted conditions, suggesting its involvement in cell division regulation under nutrient fluctuations. Conclusions: This work provides a simple and efficient method for the generation of lncRNA knock-out mutants in P. tricornutum. Furthermore, this approach can be easily adapted to other microalgae amenable to transformation, enabling the enhancement of our understanding of these important cellular regulatory transcripts in the marine environment.

PI3K/AKT signalling is a key driver of ovarian cancer progression and an important therapeutic target. Its deregulation is usually explained by genetic changes and altered signalling at the cell surface. Here, we identify an additional layer of pathway control, in which a change in the level of a cancer-associated ribosomal RNA modification alters the translation of proteins that regulate PI3K/AKT activity. SNORA81 is a high-grade serous ovarian carcinoma-associated H/ACA small nucleolar RNA that directs pseudouridylation of 28S ribosomal RNA at position U4636. Integration of tumour transcriptomes with independent SNORA81-depletion datasets identified 22 SNORA81-linked transcripts with reversed disease-associated expression, 19 of them linked to PI3K/AKT signalling. SNORA81 depletion reduced AKT phosphorylation, impaired cell growth, migration, and altered cell-cycle distribution. Molecularly, polysome profiling following SNORA81 depletion identified 524 mRNAs with altered ribosome association, 485 of which changed independently of RNA abundance and included 29 PI3K/AKT regulators. SNORA81 depletion reduced ribosome association of ITGB4, RASGRP2, IGF1 and FGF18 transcripts, and reduced ITGB4 and RASGRP2 protein abundance. IRAK1 and SPC24 acted as downstream RNA-response genes contributing to selected SNORA81-associated phenotypes. ITGB4 expression partially restored cell growth and restored downstream RNA responses after SNORA81 depletion. These findings identify rRNA modification as a translational regulator of PI3K/AKT signalling in ovarian cancer.

Intrinsically disordered regions (IDRs) mediate molecular interactions but can also promote aberrant self-association. Using donors embedded within IDR-coding sequences, exitron-type splicing can remove these sequences, giving them both protein-coding and splice-regulatory roles. Yet IDR amino-acid sequences often diverge rapidly, raising the question of whether and how their coding sequences support IDR-removing splicing across species. Using RNA-seq data from jawed vertebrates, we identified events predicted to remove most of an IDR while preserving annotated structured regions. Across species, donor-flanking nucleotide identity exceeded expectations from synonymous alternatives encoding the same amino-acid sequence. TARDBP, which encodes TDP-43, provided a model for examining how coding sequences both encode IDRs and regulate exitron-type splicing. Sequence-divergent vertebrate TARDBP terminal regions supported this splicing in human cells. In teleosts, TARDBP ohnologs (tardbpa and tardbpb) encoded compositionally similar disordered tails but differed markedly in IDR-removing exitron-type splicing. Reciprocal donor-proximal substitutions in zebrafish minigenes redirected donor choice and altered the balance between IDR-retaining and IDR-removing transcripts. In corresponding reporters, these substitutions altered protein-product balance and cytoplasmic puncta formation. These findings show that IDR-removing exitron-type splicing can recur despite sequence divergence and that donor-proximal TARDBP sequences both encode IDR segments and regulate IDR inclusion, linking RNA processing to protein behavior.

The long-range RNA-RNA interaction between the termini of a flavivirus genome is controlled by multiple local RNA secondary and tertiary structures. However, a complete picture of this essential genome cyclization process is still elusive. Herein, by using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), SHAPE and mutational profiling, in vitro RNA binding assays, and reverse genetics approaches, we demonstrated that the hairpin 2 in the capsid-coding region (cHP2) is required for efficient replication of yellow fever virus and functionally interacts with the 3'-small hairpin (sHP) during genome cyclization. Interestingly, through a combination of multidimensional chemical probing and functional characterization assays, it was discovered that the sHP forms previously unknown tertiary interactions with an immediately upstream GACG region embedded within the 3'-cyclization sequence. This crucial tertiary RNA motif accounts for a potential negative regulatory mechanism of genome cyclization and is universally conserved among the mosquito-borne-related flaviviruses, which constitute a major phylogenetic branch within the subgenus Euflavivirus, genus Orthoflavivirus. Our study highlights the intricacy of RNA-based regulation in flaviviruses.

Pitt-Hopkins syndrome (PTHS) is a severe neurodevelopmental disorder caused by haploinsufficient expression of transcription factor 4 (TCF4). Increasing expression from the remaining functional TCF4 allele could provide a mutation-independent strategy for patients with diverse pathogenic variants. We integrated genomic datasets with strand-specific RNA analysis to identify transcribed regions of the TCF4 promoter suitable for targeting by activating double-stranded RNAs (dsRNAs). Screening of overlapping promoter-targeted dsRNAs identified a hotspot for increased TCF4 expression. Activation was retained by a lead dsRNA containing extensive 2'-fluoro, 2'-O-methyl, and phosphorothioate modifications used in therapeutic oligonucleotides. Anti-AGO2 immunoprecipitation showed association with RNA transcribed from the targeted promoter region. dsRNA-mediated activation was associated with increased RNA polymerase II occupancy. These findings establish promoter-targeted dsRNAs as chemically tractable starting points for increasing TCF4 expression and support their investigation as a mutation-independent strategy for haploinsufficient disease.

Nanopore direct RNA sequencing has emerged as a powerful technology for studying RNA and modifications at single-molecule resolution. Although originally developed for polyadenylated RNAs, protocol adaptations have extended its application to non-polyadenylated RNA species, including tRNAs and rRNAs. However, the lack of open and flexible (de)multiplexing solutions remains a limitation for large-scale sequencing studies. Existing community-driven approaches enable multiplexing of polyadenylated RNAs and, more recently, tRNAs, but do not fully address the challenges of targeted sequence-specific Nanopore sequencing of non-polyadenylated RNAs. Here, we present an improved sequence-specific direct RNA sequencing protocol for ribosomal RNAs, achieving enhanced capture and longer read lengths as demonstrated for human 28S rRNA. We complement this protocol with an expanded WarpDemuX software toolbox for flexible sequence-specific (de)multiplexing. We validate the approach using human cell lines carrying knockouts of different C/D box snoRNAs involved in site-specific 2'-O-ribose methylation of 28S rRNA. Our results demonstrate efficient multiplexed detection of sequence-specific RNA molecules and associated modification changes. While established here using rRNA, we validate our approach on another panel of non-polyadenylated RNA species, providing an open and scalable framework for targeted direct RNA sequencing and RNA modification analysis

Transfer RNA (tRNA) modifications regulate protein synthesis, yet their sequence positions and stoichiometries are known only for a subset of tRNAs in model organisms. Aplysia californica is a key neurobiological model in which neuronal modified ribonucleosides have been linked to changes in animal behavior. However, Aplysia tRNA modifications have not been mapped, obscuring their roles in modulating neuronal functions. To address this knowledge gap, we leveraged hybridization pulldowns and an integrated mass spectrometry (MS) strategy combining nucleoside profiling, intact tRNA measurements, and bottom-up tRNA modification mapping to determine the full sequences of Aplysia tRNAGluCUC and tRNALysUUU. Sub-stoichiometric modifications including dihydrouridine (D), 5-methylcytidine (m5C), 5,2'-O-dimethyluridine (m5Um) and 1-methyladenosine (m1A) on tRNAGluCUC varied across tissues, suggesting tRNA modifications are regulated to accommodate translational needs. In addition to conserved thio-modifications in the anticodon loop of tRNALysUUU, we report 2-thiocytidine (s2C) at position 32, representing the first detection of s2C in eukaryotes. We also developed a generalizable method to directly detect tRNA halves by intact LC-MS, revealing endogenous tRNA cleavage patterns. Our findings illustrate the advantages of combining intact and bottom-up MS, reveal s2C as an unexpected tRNA modification in a classical neurobiological model and offer new insights into post-transcriptional regulation of animal physiology.

Angiogenin (ANG), a member of the RNase A superfamily, regulates multiple aspects of RNA metabolism, including rRNA processing and the generation of tRNA-derived small RNAs. However, its function during early mammalian embryogenesis remains unknown. Here, we identify a previously unrecognized requirement for ANG activity during early mouse embryonic development. Pharmacological inhibition of ANG during the first embryonic cell cycle did not prevent progression to the 2-cell stage but markedly impaired subsequent development to the 4-cell stage, indicating a critical requirement for ANG activity during the maternal-to-zygotic transition. Transcriptome-wide analysis revealed that ANG inhibition disrupted the normal remodeling of the embryonic transcriptome, characterized by persistent accumulation of a subset of transcripts that are normally reduced during the 2-cell transition, including many maternally deposited RNAs, together with altered expression of transcripts associated with zygotic genome activation. Small-RNA profiling further revealed widespread alterations in miRNA and tRNA-derived small-RNA populations. Notably, a large proportion of transcripts retained following ANG inhibition contained predicted target sites for miRNAs that were concomitantly reduced, suggesting a potential link between ANG-dependent small-RNA regulation and maternal transcript clearance. ANG inhibition was also associated with altered abundance and organization of components of RNA regulatory granules, including Tudor, TIAR, Dcp1a and Ago2, and increased association of selected retained transcripts Reep4 and Esrrb with Tudor-positive structures. Together, these findings identify ANG activity as an important regulator of RNA remodeling during the maternal-to-zygotic transition and suggest that ANG coordinates maternal transcript clearance, small-RNA dynamics and mRNP organization to support early mouse embryonic development.

Background: X-chromosome inactivation (XCI) evolved independently in eutherian and marsupial mammals, where it is orchestrated by the unrelated long non-coding RNAs Xist and RSX, respectively. Xist organizes a repressive nuclear compartment through multivalent RNA-protein interactions, but whether RSX exploits similar biophysical principles remains unknown. A recent paper has identified bona fide RSX interacting proteins. Results: We integrated proteome-scale RNA-protein interaction prediction, experimental validation, phase-separation propensity analysis, functional annotation and comparative RNA-structure modelling to characterize the RSX interaction landscape. Using catRAPID, we ranked 1,168 RNA-binding proteins from the native Monodelphis domestica proteome. Predictions were significantly enriched for experimentally identified RSX interactors, with 4.85-fold enrichment among the top 50 candidates (P approximately 1.6 x 10-5), increasing to approximately eightfold for proteins shared by the experimental RSX and Xist interactomes (P approximately 2 x 10-6). Among 30 high-confidence RSX interactors, 13 were experimentally supported, 17 were previously unrecognized candidates and 17 exhibited high phase-separation propensity. The network was enriched in ribonucleoprotein granules and nuclear bodies and converged on m6A regulators and SR-family splicing factors. Comparative modelling detected no conserved secondary or tertiary architecture between RSX and Xist. Conclusions: RSX and Xist appear to have converged not through RNA sequence or global structure, but through recruitment of related, condensation-prone protein networks. These findings identify interaction-network and biophysical convergence as a potential principle of lncRNA-mediated chromosome regulation and provide testable candidates for determining whether RSX establishes a condensate-like compartment on the marsupial inactive X.

Bacterial cultures entering stationary phase (SP) undergo complex physiological alterations, helping the cells to survive when medium resources are exhausted. The onset of various stress responses underlying SP physiology is generally believed to be due to reactions of individual cells to the environmental cues such as starvation, toxic metabolites, etc. The SP physiological state makes bacteria unsuitable for the replication of most bacteriophages; however, some phages are able to infect and multiply in them. We investigated Escherichia coli phage DH23 growth in SP cultures of E. coli MG1655 infected at different hours post inoculation (ages). Under our conditions the cells enter SP at about 6 h, but the phage replication was possible till 11h and then dropped abruptly by 13h of culture aging, indicating an abrupt physiological change to deeper dormancy during SP. The onset of this phage-inhibiting middle-SP dormancy (MSPD) turned out to be mediated by intercellular communication mediated by the middle-SP signal particles (MSP) which are larger than 100 kDa and contain both RNA and DNA. These MSP are inactivated by RNAse or by DNAse, enabling phage growth in 14h-old cultures and lifts the tolerance of such cultures to kanamycin. Moreover, the RNAse or DNAse treatment makes 24h culture spent media suitable for additional cell growth, and cultures initiated in fresh LB supplemented with RNAse reach an OD600 about 1.5 times higher compared to the untreated control. This indicates that MSP signaling influences cell physiology well before the MSPD onset and helps the population cease growth pre-emptively to save about 1/3 of medium resources to support the viability during SP.

RNA structure governs the function of non-coding RNAs and influences mRNA stability and translation, making testing of models of experimental RNA structure relevant to fields spanning structural biology, virology, and molecular therapeutics. Here we present FAST-MaP (Fast and Accessible Sequencing Technology for Mutational Profiling), a protocol that enables per-nucleotide RNA structure characterization using standard molecular biology equipment and requiring no sequencing infrastructure or bioinformatics expertise. RNA is chemically modified with orthogonal probes (2A3 and DMS), reverse-transcribed, and PCR-amplified to produce dsDNA amplicons that are submitted directly to a commercial primer-less sequencing service. Returned FASTQ files are processed through a freely accessible web server to generate normalized reactivity profiles within minutes. The complete protocol, from DNA template to structural data, can be completed in approximately one week. We illustrate the workflow on a 659-nucleotide RNA, demonstrating how to test structure preservation across buffers, and how to test specific secondary and tertiary structure predictions of the RNA from computational modeling or cryo-electron microscopy. The protocol requires only standard molecular-biology skills and does not require sequencing or bioinformatics expertise.

Minor introns are rare but evolutionarily conserved introns in the human genome. Defects in minor spliceosome components have been reported to cause several human diseases. Mutations in ZRSR2 are frequently observed in myelodysplastic syndrome (MDS) patients. However, the precise molecular mechanism by which these mutations cause a loss of ZRSR2 function in minor intron splicing remain elusive. In this study, we characterized the impact of MDS-causing missense mutations in ZRSR2. Most of the mutant ZRSR2 proteins we tested showed different subnuclear localization pattern from that of wild-type ZRSR2 protein. Co-immunoprecipitation (co-IP) analysis revealed that UHM missense mutations disrupted the interaction between ZRSR2 and SF3B complex. We also demonstrated that ZRSR2 preferentially immunoprecipitated minor spliceosome components over major spliceosome components and it directly interacts with SF3B1. Furthermore, co-IP analysis utilizing a variety of ZRSR2 mutants revealed that mutations in the UHM and Zn2 domains significantly attenuated the interaction with SF3B complex. Finally, UHM and Zn2 mutations markedly disrupted the supporting activity of ZRSR2 for minor intron splicing. Our results indicate that disruption of the physical interaction with SF3B complex is one of the hallmarks of MDS-causing mutations of ZRSR2, which results in aberrant splicing of minor introns.